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ABSTRACT 

The normal parameters are a non-linear transformation of the cosmological parameters whose like- 
lihood function is very well-approximated by a normal distribution. This transformation serves as an 
extreme form of data compression allowing for practically instantaneous calculation of the likelihood of 
any given model, as long as the model is in the parameter space originally considered. The compression 
makes all the information about cosmological parameter constraints from a given set of experiments 
available in a useable manner. Here we explicitly define the normal parameters that work for the current 
CMB data, and give their mean and covariance matrix which best fit the likelihood function calculated by 
the Monte Carlo Markov Chain method. Along with standard parameter estimation results, we propose 
that future CMB parameter analyses define normal parameters and quote their mean and covariance 
matrix. 

Subject headings: cosmology: theory — cosmology: observations — cosmic microwave background — 
methods: data analysis 



1. INTRODUCTION 

The challenge of turning a CMB dataset into cosmolog- 
ical parameter constraints is one that has been solved by a 
series of data compressions. The information in the time- 
ordered data is compressed into a map (Wright et al., 1996; 
Tegmark, 1997a). The information in the map is then com- 
pressed into a power spectrum (Tegmark, 1997b; Bond 
et al., 1998, 2000; Wandeh et al., 2001; Bartlett et al., 
2000; Hivon et al., 2002). Finally, the information in the 
power spectrum is then compressed into cosmological pa- 
rameters (Lineweaver et al., 1997; Benoit, 2002). 

This last step, however, is more of a data "explosion" 
than a data compression. Although the number of param- 
eters is indeed small (^ 10) the non-Gaussian distribution 
of their errors may be characterized by 10^° numbers for 
grid-based likelihood evaluation (Tegmark & Zaldarriaga, 
2000) or perhaps as little as 30, 000 for the Monte-Carlo 
Markov Chain method (Christensen et al., 2001). 

This explosion has undesirable consequences. The full 
constraining power of the data is not used by the cosmo- 
logical community. Typically authors report not the full 
(cumbersome) likelihood, but projections or marginaliza- 
tions of it down to one or at most two-dimensional spaces. 
Given the near-degeneracies that exist in the probability 
distribution, and its non-Gaussian nature, these final steps 
lose significant information. 

Our work here provides the final step of compression 
in the data analysis pipeline. We compress the probabil- 
ity distribution of the cosmological parameters down to 
the ^ 50 parameters of an analytic form for the likeli- 
hood (a mean and covariance matrix). This step makes 
the full information in the parameter likelihood function 
easily useable. 

Our procedure is analogous to the 'radical compres- 
sion' (Bond et al., 2000) used for compressing the uncer- 
tainty in the power spectrum estimates, whose distribtu- 
ion is also non-Gaussian. In both cases a non-linear vari- 
able transformation is used, with the property that the 
transformed variables are well-approximated by a normal 



distribution. 

We were inspired to search for a set of normal parame- 
ters by a proposal that the C; can have a linear dependence 
on a set of well-chosen parameters that span the space of 
possible models (Kosowsky et al., 2002). Clearly such a 
parameter set would be tremendously valuable for param- 
eter estimation since the C; for any given model could be 
calculated practically instantaneously simply by summing 
the terms in the first order Taylor expansion. Such a set 
would also, given sufficiently Gaussian errors on Ci, have 
errors with a normal distribution. 

Although inspired by Kosowsky et al. (2002), our goal 
here is entirely different. It is to find a set of parameters 
(the normal parameters) for which the likelihood is Gaus- 
sian. Although the approximate linearity demonstrated by 
Kosowsky et al. (2002) gave us confidence to try to find a 
set of normal parameters, it by no means guaranteed suc- 
cess. The linear-response approximation has not yet been 
demonstrated to be sufficiently accurate for purposes of 
parameter estimation given MAP data (and certainly not 
given current data) and the errors on Ci are not actually 
Gaussian. 

While Kosowsky et al. (2002) are developing a tool for 
parameter estimation from CMB data, our work is aimed 
at being able to easily take full advantage of such param- 
eter estimation. The most straightforward application of 
CMB data compressed in this form is for simultaneous 
analysis of CMB constraints with those from other cosmo- 
logical probes. In the past to do this analysis one would 
have to entirely reproduce the reduction from bandpowers 
(plus their window functions, offset log-normal parame- 
ters. Fisher matrices, and calibration uncertainties) to cos- 
mological parameter likelihoods. This procedure requires 
the organization of a lot of data as well as hundreds of 
thousands of angular power spectrum calculations. In con- 
trast, with the compression to normal parameters in hand 
one need only perform a variable transformation and eval- 
uate a 10-dimensional Gaussian. We use a similar set of 
parameters as KMJ, and demonstrate explicitly the valid- 
ity of the normal approximation by numerical computation 
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of the likelihood given current CMB data. 

One can also simplify a probability distribution by choos- 
ing linear combinations of the parameters that diagonal- 
ize the covariance matrix, the so-called parameter eigen- 
modcs (Efstathiou & Bond, 1999). However, the resulting 
parameters are still dependent and non-Gaussion if the 
distribution of the original parameters is non-Gaussian. 
We do find diagonalization of the covariance matrix oj the 
normal parameters to be useful, both for our fitting pro- 
cedure and for gaining an understanding of the nature of 
the constraints the data places on the parameter space. 

We test our procedure on current data and provide the 
mean and covariance matrix to the normal parameters that 
best fit the likelihood given current data. As we go to 
press, this fit is now out of date due to the data from the 
Wilkinson Microwave Anisotropy Probe ( WMAP) ^ . How- 
ever, we expect these data to further improve the validity 
of the normal approximation to the likelihood and thus 
our method will be of lasting value. 

In section II we describe the normal parameters. In sec- 
tion III we describe our calculation of the exact likelihood 
via the Monte Carlo Markov Chain method (Christensen 
et al., 2001). In section IV we describe our procedure for 
finding the mean and covariance matrix that provide the 
best fit to our MCMC-calculated likelihood. In section V 
we present our results. In section VI we discuss applica- 
tions and conclude. 

2. THE NORMAL PARAMETERS 

The normal parameters are a non-linear combination of 

the six cosmological parameters that we consider here - 
wt, wd, i^A, z, A= (2.73 X 10^)2 Mpc~^P$(0.05Mpc"^) 
and Tis (baryon density, dark matter density, dark energy 
density, reionization redshift, amplitude and scalar tilt of 
primordial power spectrum respectively). In particular, 
we replace f^A, A and by 

~ ef = rs/dA (1) 

'0.05Mpc-i\'-"'^_2, (2) 
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where is the soimd horizon at recombination, (Ia is the 
angular diameter distance to the recombination surface 
and we take fcpivot = 0.067 Mpc~^. The distinction be- 
tween 6s and Of is explained below. The numerical fac- 
tors are chosen so that A* is of order imity. We use t to 
denote the optical depth to Thomson scattering from here 
back to some time before reionization and after recombi- 
nation. The resulting variable set, {wi,,^^, 6^, ^;}, 
has a probability distribution, given CMB data, that is 
well- approximated by a Gaussian as we will demonstrate 
below. 

We chose the normal parameters by considering what 

combinations of cosmological parameters affect the fea- 
tures in the angular power spectrum over the range in 
which the data are highly constraining (100 < I < 500). 
The parameter A* describes the overall amplitude, t is 
chosen to correlate with the ratio of the amplitude of the 
second peak to the first peak (at fixed w^)) and scales 
the features horizontally. 

^ http://map.gsfc.nasa.gov 



Hu et al. (2001) have investigated this phenomenology 
as well. For G5 we use their approximation to the angular 
size of the sound horizon, which has the advantage of be- 
ing an algebraic expression^, rather than the exact sound 
horizon. Of'. This distinction is important because they 
can differ by more than the uncertainty in ©5. Hu et al. 
(2001) also define a parameter H2 analogous to our t but 
with a different parameter dependence. 

The importance of O., for understanding the Ci is widely 
recognized (Efstathiou & Bond. 1999; Tegmark & Zal- 
darriaga, 2000; Hu et al., 2001; Kaplinghat et al., 2002; 
Kosowsky et al., 2002). It is one of the best-determined 
cosmological quantities: Og = 0.°59 ± 0.°01 (Knox et al., 
2001). 

The parameter A* sets the overall amplitude of the spec- 
trum at / > 100. At these t values, Thomson scattering 
depresses the amplitude by exp(— 2r). The value of fcpivot 
was chosen to decorrelate A* and ns- In an idealized case 
with uniform relative error bars on l{l -\- l)Ci/(27r) from 

some /min to /max WC WOuld CXpCCt fcpivot = Vmin/max/?70 

where % = (14±0.6) Gpc is approximately the coordinate 
distance to the last-scattering surface. Taking /niin = 10 
and /max — 1000 we expect fcpivot ^ 0.007Mpc^^. At fixed 
T (or z) we find this value works well. Allowing r to vary 
leads to a correlation with ns- A result of this correla- 
tion is that decorrelation between A* and ng is best done 
with a much larger fcpivot of about 0.067 Mpc~^. The r- 
ns correlation arises because both rely considerably on C; 
measurements at / < 100. 

The least "normal" of the normal parameters is z, which 
is not very well-constrained by CMB data. The situation 
may improve quite soon with large-angle polarization data 
from MAP (Kaplinghat et al., 2002). Or, it is already 
not a problem if one interprets observation of a Gunn- 
Peterson trough in a 2; = 6.3 quasar as indicating z-rx ~ 6.3 
(Becker et al., 2001; Fan et al., 2002). We caution that 
though the data indicate the amount of neutral Hydrogen 
in the inter-galactic medium is rising from zero, how much 
the fraction of free electrons is increasing with z is highly 
unconstrained (Kaplinghat et al., 2002). 

We assume a step function transition of the ionization 
fraction from to 1 at redshift z. For the flat models we 
consider this results in an optical depth of 
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(Hu & White, 1997; Kaplinghat et al., 2002). We settled 
on z as a normal parameter after trial and error. We have 
found z to be more normal than either r or exp — 2r. 

We have not included tensor perturbations, curvature or 
dark energy models with w = p/p ^ — 1 in our analysis and 
so our results strictly only apply with these assumptions. 
Including all these variations (at fixed O^) would only alter 
the C; at / < 60. Their inclusion would therefore affect our 
constraints on z and ns (and therefore t), but not A*, Qg, 
u>b or u)d- 

3. LIKELIHOOD CALCULATION 

What we want to know is, given the data and any other 
assumptions we make about the world, what is the proba- 
bility distribution of the parameters? This posterior prob- 

^ We set ©s = tt/Ia where Ia is given by their equations A3-A5. 
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ability distribution can be calculated by use of Bayes' the- 
orem which states: 

P{e\d) OC P{d\0)Pprior{e) . (5) 

where d refers to data and the proportionality constant is 
chosen to ensure / P{9\d)d9 = 1. With a uniform prior 

this simply reduces to P{6\d) oc P{d\9). This probability 
of the data given the parameters is, when thought of as a 
function of the parameters, called the likelihood, £.{0). 

Often we are interested in the posterior probability dis- 
tribution for one or two parameters alone. This marginal- 
ized posterior is given by integrating over the other pa- 
rameters. For example: 

P{eue2\d) = j W}^^d6iC{6)P^,,,,{6) (6) 

where n is the number of parameters. We use the prior to 

incorporate non-CMB information such as that the red- 
shift of reionization must be greater than 6.3 (Becker et al., 
2001). 

In the following subsections we discuss first how we eval- 
uate the likelihood function at a single point and then how 
we evaluate it over a large parameter space and produce 
marginalized posterior distributions. 

3.1. Likelihood evaluation 

Here we take the data to be the measured averages of 
the angular power spectrum (called bandpowers), Df. The 
expected signal contribution to Df is given by an average 
over the power spectrum: 

Dt = Y,ul^i)fiiCi{0) (7) 

where ^^(i) is the calibration parameter for dataset a and 
Ci = l{l + l)C;/(27r) is the angular power spectrum. 

The uncertainty in the D, is non-Gaussian but well- 
approximated by the offset log-normal form of Bond et al. 
(2000). Specifically, C = exp(-xV2) where 

x' = J2 i^t - Zf) Mfj [Z] - Zf) + xU,; (8) 

Zf^\n{Dt + x,)- (10) 
Zl^MDt+Xi): (11) 
= Mij {Df + x^) {D'^ + Xj) no sum; (12) 

where M^j is the weight matrix for the band power data 
Df. The cxpt label is for experimental parameters. These 
include the calibration parameters and a beam-width pa- 
rameter, b, for the one experiment with significant and 
quantified uncertainty in their beam width. For simplic- 
ity, we take the prior probability distribution for the ex- 
perimental parameters to be normally distributed. Since 
the datasets have already been calibrated, the mean of the 
calibration parameters is at = 1. The calibration pa- 
rameter index, a, is a function of i since different power 
spectrum constraints from the same dataset all share the 
same calibration uncertainty. 

We include bandpower data from Boomerang, the De- 
gree Angular Scale Interferometer (DASI; Halverson et al., 



2002), Maxima (Lee et al, 2001)), the COsmic Background 
Explorer(COBE; Bennett et al., 1996), VSA and CBI. The 
weight matrices, band powers and window functions are 
publicly available for both VSA and DASI. For COBE we 
approximate the window functions as tophat bands; all 
other information is available in Bond et al. (2000) and in 
electronic form with the DASh package. For Boomerang, 
CBI and Maxima we approximate the window functions 
as top-hat bands, the weight matrices as diagonal and the 
log-normal offsets, x, as zero. The Boomerang team report 
the uncertainty in their beam full-width at half-maximum 
(fwhm) as 12.9 ± 1.4 minutes of arc. We follow them in 
modelling the departure from the nominal (non-Gaussian) 
beam shape as a Gaussian. For the Boomerang Df, we 
replace C; in Eq. 7 with C; exp(— ^^6^). We do not allow 
the fwhm to go below 11'. 5 or above 14'. 3; i.e. we give 
such fluctuations zero probability. To reduce our sensitiv- 
ity to beam errors, for Boomerang and Maxima we only 
use bands with maximum ^-values less than 1000. For CBI 
we use their broader mosaic'ed fields with bandpowers that 
extend to I = 1900. 

3.2. Exploring the Parameter Space 

Our first step in exploring the high-dimensional param- 
eter space is the creation of an array of parameter val- 
ues called a chain, where each element of the array, 6, 
is a location in the n-dimensional parameter space. The 
chain has the useful property once it has converged that 
P{9 e R) = N{0 e R)/N where the left hand side is the 
posterior probability that 6 is in the region R, N is the to- 
tal number of chain elements and N{9 G R) is the number 
of chain elements with 6 in the region R. Once the chain 
is generated one can then rapidly explore one-dimensional 
or two-dimensional marginalizations in either the original 
parameters, or in derived parameters, such as to. Calcu- 
lating the marginalized posterior distributions is simply a 
matter of histogramming the chain. 

The chain we generate is a Monte Carlo Markov Chain 
(MCMC) produced via the Metropolis-Hastings algorithm 
described in Christcnscn et al. (2001). The candidate- 
generating function for an initial run was a normal dis- 
tribution for each parameter. Subsequent runs used a 
multivariate-normal distribution with cross-correlations 
between cosmological parameters equal to those of the pos- 
terior as calculated from the initial run. 

The covariance matrix of the candidate-generating func- 
tion, Cg is actually a scaled version of the posterior co- 
variance matrix, Cp: Cq = aCp. Proper choice of a is 
important if the movement through parameter space is to 
be efficient. If a is too low, the acceptance rate will be 
high, but the typical step size will be small and a full sam- 
pling of the parameter space will take a long time. If a 
is too high then most steps will be to regions of greatly 
reduced probability, the acceptance rate will therefore be 
low and once again a full sampling of the parameter space 
will take a long time. As a useful measure of the speed at 
which we sample the parameter space, we define the step 
distance to be: 

Aa„,„_i = ^^A^iCpj^.A^,-; (13) 
where A6i={6n)i — {9n-i)i, n is the iteration number and 
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i enumerates the parameters. In other words we are using 
the inverse covariancc matrix as the distance metric for 
the parameter space. We automatically adjust a in the 
first 20,000 samples of a run to maximize this average step 
distance. We find it has a fairly broad plateau between 0.3 
and 0.7. In practice we calculate the square of the distance 
with the slightly simpler expression: ^ii^Oi)^ /Cp^u. 

All of our results are based on an MCMC run con- 
sisting of 754,000 iterations. For the "burn-in" the ini- 
tial 20,000 samples were discarded, and the remaining set 
was thinned by accepting every 25th iteration, resulting in 
nearly 30,000 samples. The acceptance rate was 43% with 
a = .37 and an average step distance of 0.25 indicating 
highly efficient movement through the parameter space. 
We used the CODA software (Best et al., 1995) to con- 
firm that the chain passed the Referty-Lewis convergence 
diagnostics and the Heidelberger- Welch stationarity test. 

4. FITTING TO THE LIKELIHOOD 

We do not fit directly to the ten-dimensional likelihood 

or even to its marginalization down to the six cosmologi- 
cal parameters. The reason is that in the high-dimensional 
spaces, the likelihood is very sparsely sampled. For a fairly 
coarse grid with 8 steps in each direction, there would be 
6^ > 250, 000 grid points and yet we have only 30,000 sam- 
ples. Only with the marginalizations down to even lower- 
dimensions does the likelihood become densely sampled. 

Instead we fit to the two-dimensional marginalized like- 
lihoods and from these fits reconstruct the full result. We 
typically use 50 bins in each parameter which results in 
about 100 samples in the most-likely of the 2500 bins. 

We do this fitting in the normal parameter space. That 
is, at each point in the chain we calculate the value of 
the two normal parameters of interest and histogram that 
point accordingly. For every such pair of normal parame- 
ters we find 6i and F^j^ that give the best-fit two-dimensional 
Gaussian likelihood where 

Prom these n x (n — l)/2 (= 15 for n = 6) fits we recon- 
struct the six-dimensional Gaussian. We exploit a prop- 
erty of Gaussian distributions that the inverse of Ff,-^ gives 
the elements of the covariance matrix in the full higher- 
dimensional space: Cjj = {F'^'^)~K Thus we uniquely de- 
termine the 15 independent off-diagonal elements of the 6 
by 6 covariance matrix. The diagonal elements are over- 
determined with five estimates of Cu for each i; we average 
these together for the final estimate. The full covariance 
matrix is then inverted, giving the 6-dimensional Fisher- 
matrix of the likelihood function. The 6i are also averaged 
to produce an estimate of the likelihood maximum. 

Note that each of the off-diagonal elements is obtained 
from a 2D gaussian fit to different variable pairs. There- 
fore there is no guarantee that the covariance matrix built 
from these 2x2 blocks will necessarily be positive-definite 
(which a true covariance matrix should be). We employ an 
iterative procedure to solve the above problem. We first in- 
crease the diagonal terms of C until C is positive-definite. 
We then invert C to get the Fisher matrix F = and 
find its eigenmodes. Now we identify these eigenmodes as 
the "new normal parameters" and perform the 2D fit pro- 
cedure again (this step is not computationally intensive). 



C^d oc exp 
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The resulting covariance matrix (from the 2D fits) will be 
close to diagonal and hence positive definite. We repeat 
this procedure typically about 15 times by which time the 
changes in the covariance matrix are less than 3%. We 
have also verified that the final result is insensitive to the 
arbitrary adjustments to the covariance matrix (to make 
it positive-definite) in the first step. 

The end result is an approximate likelihood, , with 
most likely values 0i and covariance matrix C^-: 

-2lnjC'' = {9,~9,)Cr}{e,-e,). (15) 

If we denote the cosmological parameters with c and the 
normal parameters with 9, the likelihood of the cosmolog- 
ical parameters is given by 



£^(^(c))lS 
o6 



where 
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(16) 



(17) 



is the Jacobian of the variable transformation. 

5. RESULTS 

In Table 1 we show the parameters of that make it a 
good approximation to our MCMC-calculated likelihood. 
For ease of interpreting the magnitude of the off-diagonal 
elements we have shown the correlation matrix which is 
given by Cij / ^CuCjj instead of Cij itself. 

We compare the fit to the exact likelihood in Fig. 1 
which shows 6 of the 15 possible marginalizations down 
to two dimensions. The marginalization of the normal fit 
is done by a Monte Carlo process. We use the normal fit 
to rapidly generate a chain whose elements are samples 
from the normal distribution. We can then manipulate 
this chain (e.g., to marginalize down to 2 dimensions) just 
like we manipulate the Markov chains. 
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Fig. 1. — Comparison of the two-dimensional marginalized pos- 
terior probability distributions of the normal parameters for the 
normal fit (shading) and the MCMC-calculated likelihood (con- 
tours). Contours and shading show where the probability is above 
exp(— 2.3/2) and exp(— 6.17/2) of the maximum value. 
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All six marginalizations down to one dimension are shown 
in Fig. 2. One can see that the fit provides an excellent 
match to the exact one-dimensional likelihoods. 

Note that although the fit is Gaussian in the six-dimensional 
space, the priors implicit in the marginalization down to 
one dimension result in slightly non-Gaussian one-dimensional 
likelihoods. This is because the marginalization is not 
done over the full parameter space but has the constraints 
2; > 6, < < 0.9, 0.45 < h < 0.95 and 0.02 < 0„ < 
1.1. 
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Fig. 2. — Comparison of the one-dimensional marginalized poste- 
rior probability distributions of the normal parameters for the nor- 
mal fit (solid, red) and the MCMC-calculated likelihood (dotted, 
black). 
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Our normal fit to the Hkelihood can also be used to 
obtain marginalized constraints on the cosmological pa- 
rameters. Marginalizations down to two dimensions and 
one dimension are shown in Figures 3 and 4 respectively. 
Note that the highly non Gaussian distributions of the 
cosmological parameters are described well by the fit that 
is Gaussian in the normal parameters. 





Fig. 4. — As in Fig. 2 but for other parameters. 
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The eigenniodcs and eigenvalues A are shown in Table 2. 
These are the eigenmodes of the fractional Fisher matrix, 
F/^ = 0iF°,6,j where Oi are the best- fit values of the nor- 
mal parameters. We look at the fractional Fisher matrix 
so its eigenvalues quantify the fractional error and not the 
absolute error. The eigenvectors are ordered from highest 
eigenvalue to lowest eigenvalue. We see that current CMB 
data constrain four parameters at the 10% level or better. 
As expected, the best-constrained eigenvector is predom- 
inantly ©s and the worst-constrained is almost entirely 
z. The second-from- worst is Ud, though even it is farily 
well-constrained with = 0.16. Eigenvectors 2, 3 and 
4 are mostly t, A* and wt, respectively, though each with 
significant amounts of other parameters mixed in. 

We have also explored the likelihood distribution of the 
parameter set proposed by Kosowsky et al. (2002). In 
2-D projections, most pairs of parameters do look Gaus- 
sian. An exception is any pair involving their variable 
Z = exp(— 2t), which is significantly less Gaussian than 
the variable we have used to parameterize reionization: the 
redshift of reionization, z. 

6. DISCUSSION AND CONCLUSIONS 

We have defined a set of parameters that, given current 
data, have a likelihood that is normally distributed. Of 
course, our simple fit to the likelihood surface will very 
soon be out-of-date due to the expected MAP data. Still, 
we expect the technique to continue to be useful. Gen- 
erally, as data become more constraining, the Gaussian 
approximation becomes better. 

The challenge in the future will be to provide good pa- 
ramcterizations for those parameters that remain poorly 
determined. Some of these only have their effects at low 
such as the parameters governing reionization. For these it 
may be useful to define parameters for which In Ci is linear, 
since the errors in this quantity are very nearly normally 
distributed at low t Bond et al. (2000). Other physical 
effects, such as gravitational lensing, are only important 
at high I where, if instrument noise is the dominant con- 
tribution to error in Ci, then C; will be nearly normally 
distributed. For these one would want to use parameteri- 
zations for which C; is linear. 

An analytic form for the likelihood makes the actual 
application of CMB parameter constraints to cosmologi- 
cal questions much easier. One application is to quickly 
calculate constraints on various combinations of the cos- 
mological parameters such as cts (Holder, 2002) or the age, 
to (Knox et al., 2001). Another is to combine the CMB re- 
sults with those from other probes to improve constraints 
or search for inconsistencies (Wang et al., 2002). Or one 
can forecast expected constraints from a combination of 
the CMB and other probes such as supernovae (Frieman 
et al., 2002) or galaxy cluster number counts (Levine et al., 
2002). Although these applications are possible without an 
analytic form for the CMB likelihood, such a form greatly 
simplifies the analysis. 

Use of normal parameters will also improve the efficiency 
of Monte Carlo evaluations of the likelihood, by providing 
an analytic generating function that is well-matched to 
the actual likelihood. Although the statistical properties 
of the post-convergence chain do not depend on the gener- 
ating function used, the time required for convergence de- 



pends critically on the generating function. In particular 
the chain rapidly converges if the generating function ap- 
proximates the likelihood well. Clearly, this issue becomes 
more important for larger parameter spaces. A practical 
way to implement the use of normal parameters in gen- 
erating Monte Carlo Markov chains would be to create a 
small chain, derive the normal parameters and the gaus- 
sian approximation to the likelihood, and then use this 
gaussian likelihood as the generating function. 

In this paper we have explicitly given an analytic form 
for the likelihood (given a certain data set) of cosmological 
parameters and shown it to be an excellent approximation 
to the MCMC-calculated likelihood. We view this as the 
final step in a long data-reduction chain that begins with 
the time-ordered data. 

We thank N. Christensen, A. Kosowsky and R. Meyer 
for useful conversations. This work was supported in part 
by NASA grant NAG5-11098. 
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Correlation matrix, rms and mean of the normal fit 

TO THE likelihood. 



eigenvector 






z 




A* 


t 


A-i/2 


1 


0.02327 


0.01554 


-0.004625 


0.9852 


-0.1152 


0.1240 


0.0090 


2 


0.3664 


0.08785 


-0.01180 


-0.1448 


-0.1890 


0.8951 


0.015 


3 


-0.2624 


-0.1551 


0.03314 


0.07245 


0.8923 


0.3232 


0.048 


4 


0.8759 


0.1070 


0.007523 


0.05664 


0.3731 


-0.2811 


0.10 


5 


-0.1706 


0.9779 


0.01899 


0.002654 


0.1190 


-0.0003358 


0.16 


6 


0.009783 


-0.01314 


0.9992 


-2.9 X 10-5 


-0.03743 


0.002542 


2.8 



Table 2 

Normal parameter. Eigenvectors and the errors on 
their amplitudes, \~^^'^ . 
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